Comparison of Different Pairing Fluctuation Approaches to BCS-BEC Crossover 



00 

o 

O 
(N 

> 

o 
O 
o 



o 
o 



iR. Levin, ^ ^Qijin Chen, ^Chih-Chun Chien and ^Yan He 

^ James Franck Institute and Department of Physics, University of Chicago, Chicago, Illinois 60637, USA and 

^Zhejiang Institute of Modern Physics and Department of Physics, 
Zhejiang University, Hangzhou, Zhejiang 310027, China 
(Dated: October 10, 2008) 

The subject of BCS - Bose Einstein condensation (EEC) crossover is particularly exciting because of its 
realization in ultracold Fermi gases and its possible relevance to high temperature superconductors. In the paper 
we review that body of theoretical work on this subject which represents a natural extension of the seminal 
papers by Leggett and by Nozieres and Schmitt-Rink (NSR). The former addressed only the ground state, now 
known as the "BCS-Leggett" wave-function and the key contributions of the latter pertain to calculations of the 
superfluid transition temperature Tc. These two papers have given rise to two main and, importantly, distinct, 
theoretical schools in the BCS-BEC crossover literature. The first of these extends the BCS-Leggett ground state 
to finite temperature and the second extends the NSR scheme away from Tc both in the superfluid and normal 
phases. It is now rather widely accepted that these extensions of NSR produce a different ground state than that 
first introduced by Leggett. This observation provides a central motivation for the present paper which seeks 
to clarify the distinctions in the two approaches. Our analysis shows how the NSR-based approach views the 
bosonic contributions more completely but it treats the fermions as "quasi-free". By contrast, the BCS-Leggett 
based approach treats the fermionic contributions more completely but it treats the bosons as "quasi-free". In 
a related fashion, the NSR based schemes approach the crossover between BCS and BEC by starting from 
the BEC limit and the BCS-Leggett based scheme approaches this crossover by starting from the BCS limit. 
Ultimately, one would like to combine these two schemes. There are, however, many difficult problems to 
surmount in any attempt to bridge the gap in the two theory classes. In this paper we review the strengths and 
weaknesses of both approaches. The flexibility of the BCS-Leggett phase and its ease of handling make it more 
widely used in T = applications, although, the NSR-based schemes are more widely used at T 7^ 0. To reach 
a full understanding, it is important in the future to invest effort in investigating in more detail the T — aspects 
of NSR-based theory and at the same time the T ^ aspects of BCS-Leggett theory. 
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I. INTRODUCTION 

The subject of BCS-Bose Einstein condensation (BEC) 
crossover has recently become an extremely active research 
area. This is due principally to the discovery |[ll |3 HI B 
H H 0, H H fl^ of superfluid phases in ultracold Fermi 
gases which exhibit this crossover. Adding to the importance 
of this work is the view espoused by a number of theorists 
ifTTl [13 [H, 13 [Hi that the high temperature superconduc- 
tors are mid-way between BCS and BEC. Now, with an un- 
ambiguous realization of this scenario in the fermionic su- 
perfluids, one has the opportunity to investigate this physi- 
cal picture more closely and, it is hoped, gain insight into 
the cuprate superconductors. Equally exciting is the oppor- 
tunity to generalize, and in the process, gain insight into what 
is arguably the paradigm for all theories in condensed matter 
physics: Bardeen Cooper Schrieffer theory. For all these rea- 
sons a large number of variants of BCS-BEC crossover theory 
have been suggested in the literature. It is the purpose of the 
present paper to present an overview of two main classes of 
theories, discussing their strengths and weaknesses. Contrast- 
ing and comparing different approaches will, hopefully, point 
to new directions for future theoretical and experimental re- 
search. 

Initial theoretical work fH,!!?!] on the subject of BCS-BEC 
crossover focussed on a ground state which was shown to be 



the same as that proposed by Bardeen Cooper and Schrief- 
fer, when it is extended to accomodate a continuous evolu- 
tion from BCS to BEC. We call this the "BCS-Leggett" state. 
Here the chemical potential is solved self consistently as the 
attractive interaction strength is varied. In this way it became 
clear that the BCS trial wavefunction was far more general 
than was originally thought. Somewhat later, Nozieres and 
Schmitt-Rink (NSR) (l^ presented a scheme for calculating 
the transition temperatures Tc, which made the case that the 
evolution from BCS to BEC was again continuous at finite 
temperature. 

The discovery of high temperature superconductivity and 
the observation that their coherence length f (or pair size) was 
anomalously small led T. D. Lee and R. Friedberg to argue that 
one should include bosonic degrees of freedom in addressiiig 
high Tc superconductors. These authors introduced [19, 231 
the "boson-fermion" model almost immediately after the dis- 
covery of cuprate superconductivity. In a similar vein, Rande- 
ria and co-workers lUlll proposed that the NSR scheme might 
be directly applicable to these exciting new materials. Subse- 
quently other theorists have applied this BCS-BEC crossover 
scenario to the high Tc cuprates dllll HHH . Additional 
support has come from the experimental condensed matter 
community among whom a number |24, 25, 26, 27] have 
presented data which can be interpreted within this picture. 
Adding to the enthusiasm is the observation of a ubiquitous 
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(albeit controversial) "pseudogap" phasej Tsi |28|] in the 
underdoped cuprates, which was argued ll29ll30ll to be consis- 
tent with a BCS-BEC crossover scenario. 

The characterization of pseudogap effects associated with 
BCS-BEC crossover was, in fact, a crucial step. It was first 
recognized that one should distinguish the pair formation tem- 
perature T* from the condensation temperature ifTll Isill . 
That the magnetic properties of the normal phase in the tem- 
perature regime between Tc and T* would be anomalous was 
pointed out on the basis of numerical calculations, on a two di- 
mensional lattice. Here it was found that the spin susceptibil- 
ity was depressed at low temperatures ll32ll and this depression 
was associated with a "spin gap" which is to be distinguished 
ll29ll from a pseudogap which affects the "charge channel" as 
well. The fact that BCS-BEC crossover theory was, indeed, 
associated with a more general form of pseudogap, thus, re- 
quired further analysis and calculations. Using the formal- 
ism of the present paper, subsequent, theoretical studies of the 
spectral function (both above I3fl] and below [33] Tc) and the 
superfluid density ll34ll showed that a normal state pairing gap 
appeared in both the spin and charge channels and, further- 
more, affected the behavior below Tc as well [34] as above. 

The BCS-BEC crossover approach was applied to the ul- 
tracold Fermi gases, by Holland and co-workers fi^ and 
by Griffin and Ohashi [36] in advance of the discovery of 
fermionic superfluidity. The two groups predicted that the 
magnetic field tuneability associated with an atomic Fesh- 
bach resonance would lead to an unambiguous realization of 
the crossover scenario. These earliest applications to ultra- 
cold Fermi gases considered a Hamiltonian rather similar to 
the "boson-fermion" model of Lee and co-workers [fl9[ I20I1 
where the bosons were related to the so-called closed chan- 
nel molecules of the Feshbach resonance and the fermions 
with the open channel. Subsequent work has shown that these 
two channel complications can be essentially ignored so that 
the description of Fermi gas superfluidity is addressed using 
the same, simpler (or one channel) model as was used in the 
cuprates. 

There is now a fairly extensive theoretical literature lfT2l[T3l 
UllH^ on the Fermi gas superfluids. Nevertheless, there are 
two main theoretical schools, which have emerged. These ad- 
dress a wide variety of different issues and experiments. The 
first of these builds more directly on the BCS-Le ggett gr ound 
state and its finite temperature extensions |(30| 39, 40ll4lll42ll . 
The second approach iHIll 37, 38, 43lHl3builds on the 
contribution of Nozieres and Schmitt Rink which addressed a 
calculation of the transition temperature Tc. The NSR scheme 
has been extended by these and other authors away from Tc 
both in the superfluid and normal phases. It is now rather 
widely accepted that these extensions of NSR produce a dif- 
ferent ground state i72l l46l [4^/ than that first introduced by 
Leggett [17] and by Eagles [Id]. This observation provides a 
central motivation for the present paper. We want to set down 
our current understanding of what is known about the NSR- 
based theories from zero to very high T and similarly, ad- 
dress how the simplest ground state of BCS-Leggett evolves 
with increasing temperature away from zero. Since the ground 
states are different, we can safely assume that the finite T be- 



havior is as well. It should be stressed that while we attribute 
these author group names to the two different schools, the 
eponymous authors are not the origin of the theory reviewed 
here. The original paper by Leggett was only concerned with 
the ground state and that by Nozieres and Schmitt-Rink reca- 
pitulated and expanded on the results of Leggett and then went 
on to compute Tc using an approach which was not associated 
with this same (BCS-Leggett) T = state. 

The task in any finite temperature crossover theory is to ar- 
rive at a characterization of the thermal excitations of both the 
normal and superfluid phases. From this analysis all transport 
and thermal properties can in principle be obtained. Without 
any detailed microscopic theory one can still anticipate the 
general features of BCS-BEC crossover theory. In the BCS 
regime and below Tc, the excitations are the usual fermionic 
quasi-particles with an excitation gap equivalent to the order 
parameter. This gap represents the energy cost of unbinding 
the condensate pairs. By contrast, above Tc this gap is absent 
and the normal state is a Fermi liquid. In the BEC regime, it 
is energetically unfavorable to break up the parrs and so the 
excitations are purely bosonic above and below Tc. In the su- 
perfluid phase, they are, moreover, gapless. In between, in the 
interesting unitary regime, the excitations are expected to be 
a mix of fermionic and bosonic character. Here, importantly 
even the normal state has some bosonic features associated 
with the formation of "pre-formed pairs". These pairs arise 
from stronger than BCS attractive interactions. As a conse- 
quence there is an excitation (pseudo)gap for fermionic ex- 
citations which appears above Tc. With progressively lower 
temperatures below Tc, more and more of these pairs drop 
into the condensate. The challenge then is to treat the strongly 
interconnected bosonic and fermionic degrees of freedom in 
the most physically correct fashion. 

Below Tc the two schools referred to above emphasize dif- 
ferent aspects of this picture. One can summarize in the sim- 
plest fashion the key differences. The NSR-based approach 
views the bosonic contributions more completely but it treats 
the fermions as "quasi-free". Many body effects are effec- 
tively absent in the fermionic dispersion relation (called E^, 
which appears in the counterpart gap equation), except via a 
renormalization of the fermionic chemical potential fj* . The 
BCS-Leggett based approach treats the fermionic contribu- 
tions more completely than the alternative approach but it 
treats the bosons as "quasi-free". Many body effects in the 
bosonic dispersion, (which we call il°), are absent in the gap 
equation except via an effective pair mass renormalization 
M*. 

More specifically, the NSR approach incorporates a lin- 
ear dispersion in the bosonic degrees of freedom at small 
wavevector which is associated with the collective mode spec- 
trum of the condensate. However, Tc is calculated in the same 
way as for non-interacting fermions, except for the renormal- 
ization in /i*. The BCS-Leggett based approach in effect ap- 
proximates the bosonic degrees of freedom associated with 
the non-condensed pairs. While the collective modes of the 
order parameter have a linear dispersion, the non-condensed 
pairs have a quadratic dispersion and represent otherwise free 
"bosons". Here, Tc is calculated in the presence of a pseu- 
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dogap so that the condensing fermionic quasiparticles have 
admixed bosonic character. 

It is reasonable to conclude that the NSR based scheme 
approaches the crossover between BCS and BEC by starting 
from the BEC limit and the BCS-Leggett based scheme ap- 
proaches this crossover by starting from the BCS limit. For 
the former, indeed, the boson-like propagators which one de- 
duces are found to have many similarities to Bogoliubov the- 
ory for true bosons. It is claimed ll48ll that the NSR-based ap- 
proach is most accurate at temperatures low compared to Tc, 
presumably because there the bosonic degrees of freedom are 
those associated with the condensate and its collective modes. 
Thus, it is reasonable to assume that this represents the better 
ground state. By contrast the BCS-Leggett based scheme is 
more suitable at moderate temperatures within the superfluid 
phase and up to the pairing onset temperature T* . Above T* , 
the two approaches can be viewed as equivalent. 

Ultimately, one would like to combine these two schemes. 
There are, however, many difficult problems to surmount in 
any attempt to bridge the gap in the two theory classes. Not 
only is it difficult to effect such a combination but, thus far, 
there is no mean field theory of a weakly interacting Bose 
gas ll49i I50I1 which addresses the entire regime from T — 
through and above Tc and which does not have at the same 
time a problematic first order transition. Thus, for the in- 
teracting Bose gases, there is no counterpart of BCS theory 
which works so well over the entire temperature range. These 
complications, in true Bose systems, appear to be transmit- 
ted to NSR based theories of the Fermi gases. These spurious 
first order transitions [5 1.1 can lead to derivative discontinu- 
ities in the density profiles at the condensate edge and non- 
monotonic or discontinuous behavior in the superfluid den- 
sity, even in the intermediate or unitary regime. BCS theory, 
by contrast, exhibits none of these effects. If one is to find a 
smooth crossover between BCS and BEC at all temperatures 
T these issues will need to be overcome. 

Additional problems appear if one tries to bridge the gap by 
starting with the BCS-Leggett based scheme. The first task is 
to establish how non-condensed pair effects modify the col- 
lective mode spectrum. This appears to be a difficult prob- 
lem. While, some progress has been made [52] towards com- 
puting pseudogap effects on the Anderson Bogoliubov mode, 
there is, however, an even greater difficulty in coupling the 
non-condensed pairs with the renormalized collective modes. 
To arrive at this hybridization, one needs to introduce boson- 
boson coupling which requires that one go beyond the sim- 
ple T-matrix scheme which one considers in addressing the 
non-condensed component. This is not to say that the cou- 
pling between condensed and non-condensed pairs is absent, 
it must be there in the ultimate theory, but it will be difficult 
to implement. 

To summarize, the ground state produced by NSR-based 
approaches is likely to represent an improvement over that of 
the BCS-Leggett based approach, particularly when it comes 
to quantitative comparisons, and most particularly when the 
system is on the BEC side of resonance. However at the 
semi-quantitative or qualitative level one is often required to 
consider the BCS-Leggett ground state, and its finite tem- 



perature extensions, since globally this state behaves more 
smoothly. Moreover, this state is easier to handle and can ac- 
comodate inhomogeneities via Bogoliubov deGennes theory 
ll53[|54ll55l ^561. It is the also the primary way to study phases 
with population imbalance L57>.58,,59,,60„,61,1 . particularly in 
the presence of a trap. 

An understanding of BCS-BEC crossover provides an ex- 
cellent vehicle for reviewing the central features of two types 
of mean field theories: strict BCS theory and the theory(s) of 
the weakly interacting Bose gas. In both systems there is the 
potential for carrying some confusion over to the crossover 
problem, since there are important "degeneracies" which are 
not general and which occur at each endpoint. In strict BCS 
theory the order parameter A^c is the same as the excita- 
tion gap A. This relationship cannot persist in BCS-BEC 
crossover. In the Bogoliubov theory of the weakly interacting 
Bose gas the collective mode frequency is the same as the sin- 
gle particle excitation energy. This degeneracy derives from 
the coupling between the order parameter and single particle 
excitation spectrum. This situation is not the case in BCS the- 
ory. The way in which the linearly dispersing order parameter 
collective modes interact with the quasi-particle (fermionic) 
excitations and the extent to which they couple is subtle in 
BCS theory. 

Indeed, if one applies the Landau criterion to a magneti- 
cally dirty, gapless superconductor, (where importantly it is 
found that A 7^ A^c) it must, of course, reveal that super- 
conductivity is stable. This Landau criterion should not, then, 
refer to all possible excitations of the system but only those 
which couple to the condensate, that is, associated with the 
density fluctuations ll62ll . The gapless single particle excita- 
tions do not compromise superfluidity and thus one can pre- 
sume that they do not couple directly to the collective modes. 
An analogous inference can, then, be made about a clean 
BCS superconductor which suggests a decoupling between 
the condensed and non-condensed components- at the strict 
BCS level. 

There is another avenue for confusion. The flexibility of 
the BCS-Leggett phase and its ease of handling make it more 
widely used in T = applications, although, the NSR-based 
schemes are more widely used at T ^ 0. One has seen just 
this dichotomy in the original paper ifTsll by Nozieres and 
Schmitt-Rink. To reach a full understanding, it is important, 
then, to invest some effort in investigating in more detail the 
T — aspects of NSR-based theory and at the same time the 
T 7^ aspects of BCS-Leggett-based theory. 

The remainder of the paper is divided into four sec- 
tions. Section inipresents a theoretical overview of BCS-BEC 
crossover theory beginning first with an alternative presenta- 
tion of strict BCS theory at general temperatures T, which 
provides general insights. Then a brief overview of the ground 
state equations for the BCS-Leggett approach is presented. 
Sections |III] and |IV] give a more detailed description of the 
BCS-Leggett and Nozieres, Schmitt-Rink theretical schools, 
respectively, at general temperatures T. There we review the 
general equations and the specific application to the BEC limit 
as well as the superfluid density. Other issues are discussed 
as well which pertain to special features of each of the two 
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schools. These sections are more technical and they can be 
skipped by a reader so inclined who is advised to go directly 
to Section [V] Section [V] summarizes crucial comparisons be- 
tween the two theoretical schools. Many of these are pre- 
sented in the form of two tables. In addition we compare plots 
of the transition temperature in a homogeneous and trapped 
configuration and plots of the density profiles. Our conclu- 
sions are summarized in Section lVll 



Note that this state represents an essentially ideal Bose gas 
treatment of the pair degrees of freedom in the sense that it 
can be written entirely in terms of a single "Bose" operator 
with net zero momentum 



^U-^^T^Lki- (3) 



II. THEORETICAL OVERVIEW 



A. Ground State Wavefunctions 



One can also contemplate something closer to a 
Bogoliubov-level wave function which can be simply written 
for the case of point bosons. A reasonable ansatz is: 



We begin with a summary of possible ground state wave- 
functions for describing BCS-BEC crossover The simplest 
one is that of BCS-Leggett 



(1) 



where and Wk are variational parameters. If we define 
Ok — Vk/u]i we may write 



*oocexp(^ak&J,fe)|0) 



(2) 



Bogoliubov/ 



^ q>0 ^ 



For the fermionic system, a natural extension, which has 
been discussed in the literature [63] can be written as 



IV'i) = exp 



K 



K 



4! ^ 

K s 



/^Ki K2 K3 K4 "^Ki "^Ka '^Ks "^Lj 



(5) 



where each represents a shorthand notation for k,(Ti, and 
— K refers to a reversal of both the momentum and spin. In 
actuality, it has been shown that to recover a consistent treat- 
ment of Lee- Yang contributions, and to include the exact con- 
straint on the inter-boson scattering length ll64ll . it is necessary 

to keep terms of the form ^ J^k-^ 7ki ■ Keclci • • • Crb- 



We stress that this Bogoliubov-based wavefunction is not 
the basis for extended NSR theories. Nevertheless, this hier- 
archy of ground states should underline the observations made 
above, that we are dealing with two different and comple- 
mentary treatments of the bosonic degrees of freedom, when 
we investigate these two different approaches to BCS-BEC 
crossover theory. It should be stressed that, despite some con- 
fusion in the literature, bosonic contributions are present in 
the BCS-Leggett scheme as well, but they are appear as less 
strongly correlated than their counterparts in the NSR scheme. 
This point is re-inforced by a discussion of the BEC limit in 
Section IIIIBI This point is also reinforced by a recognition 
of the extensive fluctuation literature in BCS superconductors 
(at low dimension), which bears strong similarity i65i1 to our 
discussion of the BCS-Leggett approach. 



B. Strict BCS Theory and BCS-Leggett Ground State 

We begin by recasting strict BCS theory in a slightly differ- 
ent way which replaces the usual Gor'kov F functions with the 
product of one dressed and one bare Green's function. This 
alternate representation builds a basis to extend to the BCS- 
Leggett phase. We define the T-matrix for a BCS superfluid 
as 

tsciQ) = -Ald{Q)/T, (6) 

where Q is a four-vector and A^c is the superfluid order pa- 
rameter. This leads to the fermionic self energy, given by 

S^^^(X)-^t,,(Q)Go(g-X) (7) 
Q 

so that T.^^^{K) = -Al^Goi-K). Here, and throughout. 
Go is the Green's function of the non-interacting system. We 
write 

g^^'^{k) = [g^\k)^j:^^'^{k)]-\ (8) 

The well known BCS gap equation is: 

l + ?7^G^^^(if)Go(-X) = 0, T<rc, (9) 

K 
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which can be written in the more familiar form 



A,,(r) = -[/5]A,e(r)^^|g^, (10) 



where U is the attractive interaction which drives superfluid- 
ity. Here 



(11) 



where ek = k"^ /2m is the bare fermion dispersion. 

Once the self energy is known, the two-body or transport 
properties are highly constrained through gauge invariance or 
Ward identities. For convenience, we work in the transverse 
gauge. The response kernel for a fictitious vector potential A 
in an isotropic system is given by 



K{Q) = --P{Q), 
m 



(12) 



where P{Q) is the current-current correlation function and we 
have J(Q) = K(Q)A. Now, following the standard proce- 
dure one uses a Ward identity to construct a consistent form 
for the correlation function P^^^{Q) = 



^E(^ + ?) [GiK)GiK + Q) 



K 



+ Al,G{K)Ga{-K)G{K + Q)Ga{-K - Q)], (13) 

where, for convenience, we have dropped the superscript 
BCS which must appear on all dressed Green's functions. 

The second term in Eq. (fTsT l is important here. One can 
represent this diagrammatically as in a "Maki-Thompson" di- 
agram. More traditionally this is written as the product of 
two Gor'kov F functions. After analytic continuation {Q 
(ri, q)) and taking 17 then q 0, this expression leads 
to the usual BCS result for the superfluid density 



,BCS 
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K 



(14) 



Another important collective feature of the BCS superfluid 
state is the dispersion VLq — cq for the Goldstone Boson which 
is given by solving 







U 



Y,[GiK)GiK + Q) 



K 



+ Al^G{K)Go{-K)G{K + Q)Goi-K - Q)].(15) 

Note that the four Green's functions in Eq. JTSl ) are very sim- 
ilar to their counterparts in the superfluid density. This under- 
lines the fact that the dynamics associated with BCS theory 
involves inter-pair interactions, but only within the conden- 
sate. 

We end this section by using this analysis to write the cen- 
tral T = equations for BCS-Leggett theory. The gap equa- 
tion is that of strict BCS theory at T=0 and the only difference 
is that it is solved in the presence of a self consistent equation 



for the fermionic chemical potential /i, which must vary as the 
attractive interaction U varies: 



A,e(0) = -C/^A,,(0)^^ 



with 



E 



EBCS 



(16) 



(17) 



Importantly, we note that an equation analogous to Eq. (fTSb 
can also be used throughout the crossover as the basis for ad- 
dressing collective behavior of the order parameter such as the 
superfluid density and condensate sound mode L5Z,i66>,i67ll . In 
the BCS regime this yields c(T = 0) = vp/^/^- , while in 
the BEC Umit c{T = 0) « {Anna b /Ml). We define the 
inter-boson scattering length qb = 2a and Mb = 2m. 

All of this is relevant to the following observations. One 
might be concerned that, since the BCS wavefunction seems 
to treat the pairs or "bosons" at a cruder level than associ- 
ated with the counterpart Bogoliubov wavefunction, that this 
quasi-ideal gas behavior would somehow destabilize super- 
fluidity. This presumption is based on the observation that an 
ideal Bose gas cannot be a superfluid. We have now seen that 
effects appearing in the collective behavior associated with the 
condensate, such as the speed of sound, do not correspond to 
those of an ideal Bose gas. We thus infer that the condensate 
can reflect a rather complex dynamics, through the effective 
incorporation of higher order Green's functions into the gen- 
eralized linear response. 



C. Characterizing the Fermionic Degrees of Freedom in 
BCS-BEC Crossover at General T 

The above summary based on strict BCS theory provides an 
underlying basis for describing the fermionic degrees of free- 
dom in both theoretical approaches to BCS-BEC crossover. 
We emphasize that the bosonic degrees of freedom are absent 
at this level and that the fermionic degrees of freedom are not 
treated in an equivalent fashion in the two theoretical schools, 
although some of the expressions representing the fermions 
look rather similar. 

We write for the "gap" and "number" equations 



^^l-2/(yO^ 



E 



mf 



mf 



(18) 



(19) 



where mf corresponds to "mean field" and the fermionic dis- 
persion is 



E;:^{T)^Jie^-^,) 



(20) 
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This system of equations has been used by both schools to find 
a reasonable estimate for the temperature at which pairing or 
the pseudogap first occurs. This is called T*, (which satisfies 
Tc < T*) and can be computed by solving for the transition 
temperature in the strict mean field equations. 

We will show that this mean field theoretic approach with 

Amf{T) = A(r) (21) 
A^iT) = Al^{T) + Al{T) (22) 

is associated with the finite temperature extension of the BCS- 
Leggett theory. We have argued that the ground state wave 
function, Eq. must necessarily also contain bosonic ex- 
citations. This can be seen most clearly when we examine 
the BEC limit in Section UlIB I Therefore, within this theo- 
retical school, one must not presume that the mean field gap 
is equivalent to the order parameter. These bosonic excita- 
tions are accomodated by decomposing into condensed 
and non-condensed contributions called A^^ and A^g respec- 
tively. We will see that the number of non-condensed pairs 
associated with the pseudogap (pg) (represented by A^^) can 
be determined once one knows their effective mass M*. And 
this, in turn, is determined by choosing a propagator for the 
non-condensed pairs for which the BEC condition on the non- 
condensed pair chemical potential, fipair ~ 0, is consistent 
with Eq. (frtl i. 

As expected, the BCS-Leggett approach, which is naturally 
associated with a T-matrix scheme, does not include all the ef- 
fects of Bogoliubov theory. Within a T-matrix approach, one 
has a choice of factoring < c^cc^c > in one of two ways: 
to yield either condensate terms A^^ or pseudogap (pg) terms 
Apg. At this level one drops terms which couple the conden- 
sate and pair excitations. To mimic the effects of Bogoliubov- 
like theory, one would need to introduce cross terms of the 
form Ag^ Apg which clearly involve higher order propagators 
and go beyond a T-matrix approach. 

By contrast the NSR-based approach uses Eq. ( fTSl l with 

A„/ = A,,. (23) 

That is, the "gap" parameter is replaced by the order pa- 
rameter. Equation ( fT9] l is not used. Rather one determines 
the fermionic chemical potential /i* by first establishing the 
bosonic propagators. The latter are taken to be the collec- 
tive mode propagators for the mf Hamiltonian but with the 
renormalized chemical potential /i*. The fermionic propaga- 
tors, which also contribute to determine /^*, are derived via a 
T-matrix approach which couples the fermions and bosons. 

From Eqs. ( 1211 1 and ( l22t we see that in the BCS-Leggett 
based approach the fermionic quasi-particle dispersion i?k, 
which appears in the gap equation, contains pseudogap ef- 
fects. That is, the fermions which pair are not the bare 
fermions. However the bosonic dispersion 0°, which also 
contributes to a separate (pseudo) gap equation, contains in- 
teraction effects in a mean field sense only via a renormalized 
effective mass M*. By contrast, the NSR-based approach is 
based on a fermionic quasi-particle dispersion in which the 
fermions which pair are the bare fermions. However, many 
body effects enter via a renormalized chemical potential /i* . 



The interacting bosonic dispersion relation fig = cq is de- 
rived. Interestingly, the complexity of both approaches, at the 
level of numerical implementation, may lie in determining ei- 
ther renormalized parameter M* or /i*, which, in a compact 
way, reflects an approximate treatment of many body effects 
in the respective theories. 



III. BCS-LEGGETT APPROACH AT FINITE T < Tc 
A. Theoretical Framework 

At issue then is the incorporation of bosonic degrees of 
freedom into the gap and number equations. The two differ- 
ent approaches build on the fact that there are two different 
ways of arriving at soft bosonic modes within a generalized 
BCS structure. These modes may come from the collective 
phase mode of the order parameter (Goldstone boson) which 
is necessarily gapless in the superfluid phase. They may also 
arise from the condition that the non-condensed pair excita- 
tion spectrum is gapless. Both of these are simultaneously 
satisfied in both theory classes. 

We turn first to the BCS-Leggett based theory, which pro- 
vides a very natural and straightforward extension of BCS the- 
ory. We note that strict BCS theory has two distinct conditions 
for soft modes of two particle propagators, one coming from 
the Goldstone boson and the other from Eq. (|9]). This observa- 
tion plays an important role in the extension of BCS-Leggett 
theory to finite T. We begin by presenting the central equa- 
tions, rather than giving a full derivation. Two of these equa- 
tions have already been written down in Section III CI for the 
superfluid regime. These are Eqs. ( fTST i and Eq. ( fT9] l, impor- 
tantly, with the substitution A„i/(T) = A(T), as in Eq. dlTT i. 

In order to quantify the pair fluctuations, our task is to de- 
compose A2(T) into AIXt) and Alg{T). The difference 
between the gap A and the order parameter Age is to be as- 
sociated with pair fluctuations (involving Apg), as should be 
implicitly evident in Eq. ( l22l i. The physical arguments which 
we apply next are rather analogous to Bose Einstein conden- 
sation: once we know the propagator for the non-condensed 
pairs we determine the number of such pairs and in this way 
we determine A^^. We can essentially anticipate the answer 
simply by counting all non-condensed pairs as 

Alg{T) = Z-^Y.^{ai,Tl (24) 

where Z is an overall coefficient of proportionality, to be de- 
termined below and 6 (w, T) = l/[exp(aj/r) — 1] is the Bose 
function. Here is the non-condensed pair dispersion. Then 
just as in BEC theory, knowing the non-condensed pair con- 
tribution (Apg) and the total (A^) one can find the condensate 
term A^^. 

To make progress we need to evaluate fi" (and Z). We 
equate the condition that the propagator for non-condensed 
pairs has zero chemical potential 

Impair — (25) 
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at and below Tc, with the gap equation Eq. ( fTSl ), where 
Eq. ( I2T] ) must be imposed, so that we have 



1 - 2/(^k) 



= 0, 



T<T, 



with 



(26) 



(27) 



Note that it is the excitation gap and not the order parameter 
which appears here. That the BCS form for the gap equation 
is equivalent to the gapless condition on non-condensed pairs 
imposes a constraint on the non condensed pair propagator 
which must be of the form 

tpg{Q) = U/[l + Ux{Q)l (28) 
where, importantly, one must take the pair susceptibility 



x(g) = 5]Go(Q-if)G(if). 



(29) 



with 



a2 =-V 

P9 



(35) 



This equation will be shown below to be equivalent to 
Eq. (EHi. We write 



(k,c.) 



(36) 



K 



from which one finds S(k, oj) w A^/[w + eu — ij] where we 
have used Eq. ( l22b . In this way one derives Eq. ( |26] |. 

Note that in the normal state (where /ipair is nonzero), 
Eq. ( l34l i is no longer a good approximation, although a nat- 
ural extension can be readily written down ll69tl . 

At small four-vector Q, we may expand the inverse of tpg 
after analytical continuation. Because we are interested in the 
moderate and strong coupling cases, where the contribution 
of the quadratic term in il term is small, we drop this term 
and thus find the following expression, which yields il° — 

/{2M*) via the expansion 



Here G and Go are the full and bare Green's functions re- 
spectively. We have met the combination GGq in the con- 
text of our review of conventional BCS theory. To expand on 
this point, note that the full Green's function is determined in 
terms of the usual BCS-like form for the self energy 



(30) 



Using this self energy, one determines G and thereby can 



evaluate tpg. The gap equation in Eq. 
tpg{Q) — 00. Similarly, using 



thus requires that 



(31) 



one derives 



E 



k 



£k - M 
Sk 



2^/(£;k) 



(32) 



which is the natural generalization of Eq. (fT9] l. 

The final set of equations which must be solved is rather 
simple and given by Eq. ( |32] |. Eq. (l24l) . and Eq. ( l26l ). This set 
has a more detailed derivation, and we summarize it by noting 
that there there are two contributions to the full T-matrix t = 
tnn + tf;r whcre tsriQ) = —^w^S(Q). Similarly, we have 



"pg + tsc where tsdQ) = — y^5{Q) 
for the fermion self energy S(-ft') = HsdK) 
J2q t{Q)Go{Q - if). It follows then that 



^pgiK) 



i;sc(k,Lj) = 



AS 

k,sc 



(33) 



A vanishing chemical potential means that tn niQ) d iverges 
at Q when T < Tc. Thus, we approximate Ollesll T.{K) 
to yield 



^pg{K)^-G^{~K)Al^ T<Tc, 



tpgiQ) 



z- 



(37) 



where Z is a residue given by 



^^pg 



on 

1 

2A2 



a=o,g=o 



2 2^/(ek 

k 



m) 



(38) 



Further details are presented in AppendixlAl 

Below Tc the imaginary contribution in Eq. ([37]) Tq 
faster than as q 0. It should be stressed that this ap- 
proach yields the ground state equations and that it represents 
a physically meaningful extension of this ground state to finite 
T. 

We note that the approximation in Eq. ( [34l i is not central to 
the physics, but it does greatly simplify the numerical analy- 
sis. One can see that correlations which do not involve pairing, 
such as Hartree terms are not included here. This is what is 
required to arrive at the BCS-Leggett ground state. It should 
be clear that, in principle, the T-matrix approach discussed 
here is more general and that in order to address experiments 
at a more quantitative level it will be necessary to go beyond 
Eq. ( [34I 1. Indeed, the simplest phenomenological correction is 
to write 



Spg(k,w) 



A^ 

k,pg 



LU + e^ — fi + ij 



Eo(k,c^). (39) 



Here the broadening 7 7^ and "incoherent" background con- 
tribution Eo reflect the fact that noncondensed pairs do not 
lead to true off-diagonal long-range order By contrast S^c 
is associated with long-lived condensed Cooper pairs, and as 
shown in Eq. ( [33] l, it is similar to Epg but without the broad- 
(34) ening. It is important to note that this same analysis has been 
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applied to describing the spectral function in the pseudogap 
iFtoI ItTIi and the superfluid phases |72] of the high tempera- 
ture superconductors, where here So(k, w) is taken to be an 
imaginary constant. 

In summary, the simplifying approximation in Eq. (|34] | 
is most problematic when the pairing gap is small so that 
other correlations and contributions (which are otherwise in 
the "background") become important. Perhaps the most noba- 
ble example of when this simplification affects the qualitative 
physics is in the population imbalanced gases. At a quanti- 
tative level, a clear shortcoming comes from the neglect of 
Hartree interaction effects. These issues are discussed in Sec- 
tion |VD] 

Finally, we present results for the thermodynamic al poten- 
tial , which is given by 

ft ^ rtf + iib, 



ilb = ^Tln(l-e-"?/^). 



(40) 



This thermodynamical potential can be used to generate the 
self consistent equations presented above 



dA 



= 



(41) 



which is equivalent to the gap equation of Eq. ( l26b Similarly, 
we have 







(42) 



which leads to the equation for the pseudogap given by 
Eq. ( |35] ). Finally, the number equation 



on 

9/1 



(43) 



which yields Eq. ( |32] |. 

We recapitulate by rewriting the central gapless condition 
for the non-condensed pairs as 



U 



1 + UJ:uGo{-K)G{K) 



oo. (44) 



tpgiO) 



This equation is equivalent to Eq. ( |26] | or Eq. (1251 1. Expanding 
tpg{Q) determines the excited pair dispersion 



n° ^ q^/2M* 



(45) 



B. BCS-Leggett Approach to EEC 



There has been some confusion voiced about whether the 
BCS-Leggett ground state requires that one ignore bosonic 
degrees of freedom. To respond (in the negative) to this con- 
cern it is useful to address the extreme BEC limit. We begin 
by making the important observation [73] that for T < Tc, 



the fermionic parameters associated with the wavefunction of 
Eq. ([T]l, A(T) and ^{T) are temperature independent in the 
BEC, for all T < Tc- This is consistent with the physical pic- 
ture of well established, pre-formed pairs in the BEC limit, so 
that the fermionic energy scales are unaffected by T below Tc- 
We now extend these qualitative observations to a more 
quantitative level. The self consistent equations in the BEC 
limit for general temperature T can then be written as 



E 

k 

E 



1 1 

2ek 2E]^ 



1 



T < Tr. 



(46) 
(47) 



where we have now introduced the usual s-wave scattering 
length, a, which is needed to regularize the gap equation for 
a contact potential. Note that we have used the T = con- 
ditions 1 17] in Eqs. ( |46] l and ( |47] i, since the Fermi function 
/(i?k) is essentially zero in the BEC limit, where E\^/T ^ 1. 
Equations ( |46] | and (|47] | are central to the BEC-theory. They 
show that even in the strong attraction limit, where the sys- 
tem can be viewed as consisting of "bosons", the underlying 
fermionic constraints on A and fx must be respected. 

It follows from the above equations that for general T < Tc, 



where the coefficient of proportionality 



(48) 



(49) 



This coefficient Z was obtained directly from the ground state 
equations ll43l[74ll . However, it can also be readily derived at 
non-zero T using the propagator for non-condensed pairs fol- 
lowing Eq. ( [38] l. Here one drops the last term involving the 
summation over free fermion states, which are clearly negli- 
gible in the BEC. That the same answer is obtained from the 
ground state and from tpg{Q) demonstrates an internal con- 
sistency of the calculations. 

We arrive at an important physical interpretation of the BEC 
limit. Even though A or Upairs is a constant in T, this con- 
stant must be the sum of two temperature dependent terms. 
Indeed it follows from Eq. (l22T i that, just as in the usual theory 
of BEC these two contributions correspond to condensed and 
non-condensed components 



condensed 
pairs 



(T) 



rioncondensed 
pairs 



(T). 



Note also that at 



noncondensed / 



(50) 



(51) 



We now rewrite the central equations ( |46] l, (|47] l in the BEC 
limit to compare more directly with the case of a weakly in- 
teracting Bose gas. 



= A^ 



(52) 
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which, in conjunction with the expansion of Eq. ( l46b . 



2m\ 
If) 



3/2 



Stt 



16 7*^ 



yields 



(53) 



(54) 



These expressions are used to ehminate the fermionic pa- 
rameters altogether and arrive at an expression which, at 
T = some have interpreted ll43l to be equivalent to 
the results of Gross Pitaevski (GP) theory. Here one iden- 
tifies an effective inter-pair scattering length qb = 2a with 
ub = n/2 which represents the number density of pairs, and 
finally hb = 2fi + h^/ ma?' is the "bare" chemical potential of 
the pairs, with Mb = 2m the pair mass. We emphasize that 
the value of 2 for the scattering length ratio is entirely dictated 
by the assumed form for the ground state, Eq. ([T]i. 

With these definitions. 



^■naBh? , . 

f^B = — TT [n). 

lUpt 



(55) 



For true bosonic systems, this GP equation is usually consid- 
ered only at T = 0, where all the pairs are condensed. In this 
regard we should interpret /i^ as a "bare" chemical potential 
which includes only a mean field Hartree shift. This is to be 
contrasted with fipair which is the chemical potential of the 
non-condensed pairs and reflects many body physics beyond 
Hartree terms. Similarly M* is the effective mass of the non- 
condensed pairs which is generally distinct from Mb ■ 

Note, however, that our derivation of Eq. jSSl l should, in 
principle, apply to all T < Tc, and, thus, the physics is 
very different from that of GP theory. Clearly /is as defined 
above is a constant in temperature. The quantity n appear- 
ing in Eq. (|55] l is, of course, temperature independent, but 
we note here that via Eq. ( fSOb it contains both condensed and 
non-condensed pairs. Their relative contribution can be deter- 
mined via an ideal gas dispersion relation with renormalized 
effective mass. This n° cx dispersion is, in turn, a conse- 
quence of the underlying gap equation Eq. ( |46] |. We stress that 
this gap equation has no counterpart in the GP theory for true 
bosons, although it can be interpreted in the fermionic context 
as reflecting the condition that fipan- — 0. 

Another essential distinction between the fermionic BEC 
and that of true bosons is that the effective mass contains in- 
teraction effects due to compositeness. The general expres- 
sion for the (non-condensed) pair mass 1/M* in the near BEC 
limit is given by 



1 



1 



ZA2 



E 



1 



21,2 



3to2A2 



(56) 



where we have used Eq. dJTl i. as well as Eq. ( |46] | and ( |47] ). 
After expanding to lowest order in na^, 



M* 



2m 1 



(57) 



Physically this increase in effective mass away from the ideal 
gas asymptote reflects the fact that pairs are less mobile, as a 
consequence of the inter-pair repulsion. This means that the 
asymptotic limit of Tc is approached/rom below, which is dif- 
ferent from the behavior found in the NSR approach 1 18]. The 
issue of whether the asymptotic limit for Tc in a mean field 
composite BEC should be approached from above or below 
has been addressed |75] in the literature, where it was argued 
in favor of the latter alternative. 

We turn now to a quantitative calculation of Tc, based on n° 
[via Eq. (|5T|l1. Equation ( ISTl ) reflects the fact that, in the near- 
BEC limit, and at Tc, all fermions are constituents of uncon- 
densed pairs. It, then, foUows that {M*Tc)^''^ (X n — const. 
which, in conjunction with Eq. $5% implies 



Tc - r° 



(58) 



Here T^ is the transition temperature of the ideal Rose gas 
with Mb — 2m. This downward shift of Tc follows the ef- 
fective mass renormalization, much as expected in a Hartree 
treatment of GP theory at Tc. Here, however, in contrast to 
GP theory for a homogeneous system with a contact potential 
ll76ll . there is a non-vanishing renormalization of the effective 
mass. 



Bogoliubov de Gennes Theory and Critical Velocity 
Calculations 



The most widely used theoretical formalism for the trapped 
Bose gases is probably Gross Pitaevski theory Jt^. This is 
because it has the flexibility to address inhomogeneous sys- 
tems and general perturbations. For the trapped Fermi gases, 
the emerging counterpart formalism appears to be Bogoliubov 
de Gennes (BdG) theory. Both BdG and GP theory are pre- 
sumed to be appropriate to the ground state. Moreover the 
ground state in question for the Fermi gases is associated with 
the BCS-Leggett wavefunction. 

The BdG equation is 



H{v) 
A(r) 



A(r) 



— E„ 



Mn(r) 
Vn{r) 



(59) 



Here H{r) = — ^. The solution of these equations is 

subject to the self consistent gap and number equations 

A(r) = -C/^u„(r)<(r)[l - 2f{E^)] (60) 



and 



n(r) = {\u^{v)\'f{E^) + \v^ir)f[l - 2f{E^)]} . 

(T,n 

(61) 

Finally, the mass current is 

J(r) = 2{Ju J{En) - J„,.[l - f{En)]}, (62) 

where J„„ = Im(u* Vu„) and J^,^ = Im(u* Vt;„). The gen- 
eral solution to the BdG equation depends on the geometries 
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and coupling constant and therefore usually requires full nu- 
merical calculation. 

This system of equations has been applied to the problem of 
BCS-BEC crossover in a number of important ways at T = 0. 
It was shown in Reference [771 that in the deep BEC, this 
scheme becomes equivalent to Gross Pitaevski theory. This 
observation may not be, in some sense, entirely surprising 
based on the arguments we have just presented in Section 
IIIIBI Moreover, one can see that there is a close analogy be- 
tween the wavefunction of Eq. ^ and that of Gross Pitaevski 
theory for point bosons. BdG theory has been used to ad- 
dress the behavior of a single vortex ISSllsilfssIl as the system 
evolves from BCS to BEC. One of the key observations here 
is that the core size (related to the coherence length i^) is non- 
monotonic with scattering length, exhibiting a minimum near 
unitarity. Moreover, there have been systematic studies based 
on BdG theory in the presence of ll56ll population imbalance. 
Finally, we want to call attention to work which addresses the 
behavior of the critical current as extracted from both vortex 
calculations ifssll and from Josephson junction studies ItsIi . 
Direct calculations using Eq. ( |62] i show a maximum in this 
current as a function of distance from the vortex core center 
and this maximum can be loosely associated with the critical 
current, Ic- Because Ic scales inversely with ^, one can in- 
fer from BdG calculations of ^ that the critical current is 
largest close to unitarity, as observed experimentally ItqIi . 

Physically, this maximum in Ic has been interpreted 
Isoll as suggesting that on the BCS side of resonance Ic is de- 
termined by the breaking of condensate pairs, while on the 
BEC side of resonance, Ic reflects the collective modes of the 
condensate. These two different mechanisms have different 
dependences on the fermionic scattering length, leading to a 
maximum which one might argue is close to unitarity. We em- 
phasize here that Ic is a property of the condensed pairs within 
BCS-Leggett theory. As noted earlier, one has to exercise cau- 
tion in applying the so-called Landau criterion in calculating 
Ic- Only those excitations which couple to the condensate 
(that is, to the density) are to be included in establishing the 
stability of the superfluid. 



D. Superfluid Density and Collective Mode Calculations 



We noted in Section IllBI that the superfluid density Us is 
highly constrained by a Ward identity once the self en ergy is 
chosen. These considerations have been applied ||2TI l39ll to 
the BCS-Leggett-based formalism where it has been shown 
that the contribution of non-condensed pairs does not directly 
contribute to a Meissner effect, as expected. The Aslamazov- 
Larkin and Maki Thompson diagrams associated with these 
finite momentum pairs cancel out and one is left with only a 
condensate contribution of the form 



\ m / V m / 



BCS 



(63) 



where [ns/m)^'-^^ is defined in Eq. (fl4] i. but with A^^ now 
replaced by A. Obviously, [ns/m)^'-^^ does not vanish at Tc, 



but because of the prefactor, the superfluid density reflects the 
order parameter and will be zero in the normal state. One can 
interpret this expression using A2^(r) = /\'^{T) ~ A^giT). 
and noting that there are two forms of condensate excitation 
which lead to a decrease in superfluid density with increas- 
ing T; the fermionic excitations, which are important to the 
extent that A(T) contains an appreciable temperature depen- 
dence below Tc, and the non-condensed pairs which enter via 

Kin 

In a related fashion, there is an extensive literature p2l, W, 
l67ll which has addressed the T = collective modes of the 
BCS-Leggett state. In the BCS limit the sound mode velocity 
is c{T = 0) = v f/VS, while in the BEC limit c{T = 0) w 
■\/ {ATmag/Mg), with the inter-boson scattering length as = 
2a, as derived in Section ITlI CI As noted in Section HTbI the 
inter-boson interactions arise in the condensate dynamics just 
as in Eq. il3[ through the presence of four Green's functions 
in the second term in this expression. 

With the introduction of non-zero temperature, the collec- 
tive mode spectrum must be deduced on the basis of a gauge 
invariant formulation of the response of the system to a ficti- 
tious vector potential, which enforces the Ward identity con- 
straints deriving from the self energy. Because A(T) ^ 
Asc{T), this calculation is much more difficult to implement. 
A lowest order approximation was discussed in Ref. lf52ll . In 
this case c(T) becomes complex, but both real and imaginary 
contributions are seen to vanish at Tc- If there is to be an even- 
tual reconciliation between the two approaches to BCS-BEC 
crossover it will be necessary, at the least to find a full solution 
to this problem. 



IV. NOZIERES SCHMITT-RINK THEORY: 
BOGOLIUBOV-BASED APPROACH TO FINITE T <Tc 

Although the normal state is similar to that originally pro- 
posed by Nozieres and Schmitt-Rink, the philosophy underly- 
ing this theoretical scheme for describing BCS-BEC crossover 
1181.1 begins with Galitskii's approach 1,821 to the dilute Fermi 
gas with repulsive interactions. Here a self-energy based on 
a particle-particle ladder is introduced. Moreover, it is clear 
that this scheme can be readily extended to the case of a weak 
attractive interaction in the normal phase, and then further 
extrapolated to the BEC limit (still remaining in the normal 
phase), where the particle-particle ladder acquires the form of 
the propagator for non-interacting bosons. It then becomes 
natural to extend this scheme to the superfluid phase, for 
which the particle-particle ladder acquires a matrix structure 
that maps onto the bosonic normal and anomalous propaga- 
tors within Bogoliubov theory. For these reasons the main 
physical emphasis was on the self-energy itself, and as a con- 
sequence on the related dynamical quantities. 

One of the virtues of this type of diagrammatic approach is 
that it is "modular" in nature, in the sense that it can be pro- 
gressively improved by including additional self-energy cor- 
rections which are thought to be important, particularly at the 
BCS and BEC endpoints. In this way, upon successive im- 
provements one can address the Popov theory for composite 
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bosons, the Gorkov and Melik-Melik-Barkudarov corrections 
fsliL etc. Of course, practical implementation of these theo- 
retical improvements suffers by the increased numerical com- 
plexity. 

For want of a better short name, we will refer to this as 
the "NSR-based approach". One can also think of it as a dia- 
grammatic T-matrix scheme which involves a matrix form of 
the T-matrix. By contrast the BCS-Leggett approach is a dia- 
grammatic T-matrix scheme which involves a scalar form for 
the T-matrix. 

In this alternative approach, Eq. (fTSl l is used iIjtI [ssll to 
yield 



l + U 



k 



1 - ViK) 

2E^ 



0, 



T <Tr 



with 



(64) 



(65) 



We can rewrite this gap equation, along with the number 
equation as 



K 

n = 2^Gii(i^). 



(66) 
(67) 



K 



Here G° is the bare matrix Green function with components 
given by = -(^ + i^n)l\{E°^? + and G?2 = 
A^c/il-E^)^ + '^n] with = fc2/2m - [i* . Note that there 
are two different levels of Green's functions which appear in 
these equations. In effect, fluctuations associated with the col- 
lective modes will appear in the number equation, but not the 
gap equation. 

The fully dressed Green's functions which include collec- 
tive mode effects are determined in terms of the matrix self 
energies 



Eii(k,w„) = -E22(-k, 

";rn(Q)G?i(Q-if), 



Q 



Ei2(k,w„) = i;2i(k,cj„) = - As 



(68) 



(69) 



from which the important dressed Green's function (which re- 
flects) the pair fluctuations and which is used in the number 
equation Eq. (l67T i can be derived [SIJ : 



1 



an (if) = Sii(if) 



G-^\K)-o^^{Ky 

El2(X)E2l(X) 



(70) 



Go^(if)-E22(if) 



HereG^i(if) = (iw„ - Ck). 

The pair propagator (which is the analogue of ipg in the 
BCS-Leggett theory) is related to the "bare collective modes". 
In particular, 



rii(Q) 



x;i(Q)x?i(-Q)-[x;2(Q)]^ 



Y^Gl^{K + Q)G\^{~K)-- (71) 

K 
K 



xUQ) = ydUK + Q)G"2ii-K), ill) 



where throughout we use the four vector notation K = 
(k, w„) and Q = {q,ilv). 

We end by recapitulating the central gapless condition of 
this class of theories: 



rii(o) = 



X?i(0) 



(73) 



The speed of sound is obtained from the finite Q general- 
ization of the denominator in Eq. (l73l i: Xii(Q)Xii(~'5) ~ 
[XuiQ)]'^ = which can be seen |80] to yield an answer 
equivalent to that obtained from Eq. ( fTSI ). Quite generally at 
small wave-vector this bosonic dispersion is given by 



= c{T)q. 



(74) 



A. Nozieres Schmitt-Rink-Based Theory in the BEC Umit: 

T<Tc 

It has been shown that fsS*] in the BEC limit, the equations 
for the collective mode propagators Fn and F12 are very simi- 
lar to the diagonal and off-diagonal bosonic Green's functions 
at the level of Bogoliubov theory |49]. In the deep BEC these 
bosonic Green's functions have a pole at 



fir 
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\2Me 



/4 



(75) 



which represents the characteristic dispersion relation for 
bosons in a weakly interacting Bose gas. Here /is is defined 
in Eq ( fSST i and AI b — 2to is the boson mass. 

The associated fermionic Green's functions are in some 
sense the more important, since these are fundamentally 
fermionic gases. In the BEC limit the equation 

Sn(k,c^„) = -^rn(Q)G?i(Q-X) (76) 
Q 

can be approximated by ignoring terms which involve Asc 
compared to |. The fermion Green's functions in the BEC 
limit are approximated as the following expressions which are 
derived in Ref. ifssll and we summarize the derivation in Ap- 
pendixlB] 



Gil = -(^k + ^^n)|[i^l + Ck + ^a 
Here the approximation 



(77) 



(78) 



has been used. This approximation is similar in spirit to that 
shown in Eq. ( [34l ). It is also demonstrated |83] that in the deep 
BEC regime, the fermion Green's function leads to 



"0 



(79) 



K 
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Here Uq and n' denote densities of condensed and noncon- 
densed pairs, respectively. Similarly, it follows that 

Gi2 = A,,/ [ujI + ^l + 2Alg + Al] . (80) 

The modified gap equation, 

A.c = C/^Gi2(if), (81) 

K 

gives fss"), in the deep BEC limit, 

^lB~(^^^^ino + 2n'). (82) 

Under these approximations, pairs in the deep BEC limit be- 
have like bosons in the Popov approximation. Although there 
is an asymmetry between the denominators of these two com- 
ponent Green's functions, one can see that the diagonal term 
has a strong similarity to the previous approach of Section 
nil A I The effective excitation gap is given by the contribution 
from condensed and excited pairs. 

To go beyond this scheme, it is necessary to incorporate 
corrections to the gap equation Eq. (|64] | with concomitantly 
those to the collective mode spectrum, so that collective mode 
effects have to be treated at a level beyond the bare modes 
of BCS theory. Some progress has been made ifssll in imple- 
menting this scheme in the BEC limit. 



B. The Controversy Surrounding the Number Equation in 
The NSR Approach 

In the original NSR approach the number equation was de- 
termined from a thermodynamical potential. Here, above and 
below 1 37] Tc, one approximates the thermodynamical poten- 
tial 

^NSR = ^mf + ^pf, (83) 

A2 



k 

- 2r^ln(l + e-^£/^), (84) 

k 

= Y.H^niQ)xU-Q)-[xUQ)n (85) 



Note that in this RPA-like scheme, only the bare pair suscep- 
tibilities x'^ are included. 

Quite generally, the number equation is given by 



(86) 



which is necessarily equivalent to Eq. $67\ . providing one 
has a complete theory. However, because of the lack of full 
self consistency, the original NSR approach was criticized by 
Serene |84]. By approximating the pair fluctuation contribu- 
tions, it corresponds to a T-matrix theory in which one takes 



only the lowest order terms in a Dyson expansion, rather than 
a full resummation, so that 



GiK) « Go{K) + Go{K)^o{K)Go{K). 



(87) 



This criticism, not withstanding, it has recently been argued 
[42,^,8^ that, below Tc one should write the number equa- 
tion as 



n ■ 



NSR. 



^ dflNSR j_ d^NSR dAsc 



dAsc dfi* 



(88) 



with dAsc/d^* determined from the BCS gap equation. 

In earlier work |37], the second term on the right hand side 
of EqdSSll was dropped. Indeed the fact that 



NSR 



dAsc 



7^0 



(89) 



in our view reflects a problem in the theory- that the gap 
equation is non-variational, or non-self consistent. This non- 
variational behavior implies that a Landau Ginsburg like anal- 
ysis, and even its generalization to first order phase transitions, 
is not possible. This anomalous term appears discontinuously 
below Tc and, it will enhance first order discontinuities at Tc, 
which may already be present in BogoUubov or Popov level 
approaches. 

Nevertheless, it has been argued that in the BEC this non- 
variational term provides a quantitative improvement over 
previous work since it evidently yields the nearly correct 18^ 
relationship between the inter-boson (as) and inter-fermion 
(a) scattering lengths. Exact few body calculations show 
that this ratio should be 0.6. It appears difficult to understand 
physically how an evidently non-self consistent gap equation 
can capture the same physics as these precise four fermion 
calculations. Indeed, this claim appears to be at odds with de- 
tailed calculations presented elsewhere which show that to ar- 
rive at this correct ratio, one must go beyond [83..87.] T-matrix 
based schemes. 

For ease in identification of these two different versions of 
NSR theory, we now refer to that based on Eq. ( |67l ) as NSR-1 
and fliat based on Eq. as NSR-2. 



C. Superfluid Density and Collective Mode Calculations 

The superfluid density Ug as a function of temperature has 
been calculated using both NSR-1 and NSR-2 like theories. 
For the former, a diagrammatic calculation of the current- 
current correlation function based on Aslamazov-Larkin and 
Maki-Thompson contributions was adopted |48], which is, in 
many ways, similar to that discussed in Section 1111 D I within 
the BCS-Leggett framework IHIll]. For NSR-2 like theories 
a framework based on changes in the thermodynamic poten- 
tial associated with a "phase twist" was adopted |88, 89]. The 
results appear to be rather similar, at a qualitative level. For 
some parameter regimes, there are either first order transitions 
at Tc or multivalued results for Ug which presumably reflect 
the analogous behavior found in Bogoliubov or Popov level 
treatments of true Bose systems 1491 pOll . See AppendixICl 



13 



An important check on these calculations is to verify that 
there is no Meissner effect in the normal state. We can follow 
the same analysis as used in Eq. (fT2l i. Quite generally, above 
Tc one has 

(-) -P:..(0) = 0. (90) 
We show below that the appropriate form for NSR-1 is 

(-) - ^Y.lS^G{K) (91) 

(a, P = x,y, z) and that this is consistent with the absence of 
a normal state Meissner effect. 

Here the current-current correlation function is 

P^AQ) = / dTe'"'"(j,,(T,q)i,(0,-q)) (92) 

K 

where A and A denote the bare and full vertices and they nec- 
essarily satisfy a Ward identity. Importantly, as shown in Ap- 
pendix |D] the two contributions cancel each other as a con- 
sequence of a Ward identity. This necessary cancellation im- 
poses an important consistency. We have presumed that the 
number equation appears as in Eq. ( |67] | (which we call NSR- 
1) which is then consistent with Eq. ( |9T| i. If, on the other- 
hand, we had assumed the number equation as in NSR-2 , 
the cancellation can be enforced as well, but only by proper 
imposition of the corresponding Ward identity. This may ex- 
plain why the results for the superfluid density in Refs. llssll 
and 189)1 were not precisely the same as those found in Ref- 
erence ll48ll . This analysis also serves to help establish those 
diagrams which must be used below Tc in order to be assured 
that there are no contributions to the Meissner current from 
non-condensed pairs. In view of the above arguments and Ap- 
pendix|D] the diagrammatic choice in Reference seems to 
be validated, although it is of interest to reformulate these cal- 
culations by explicitly imposing the Ward Identity. 

The collective mode spectrum appropriate to the NSR 
scheme was originally discussed by Griffin and collaborators 
ll37[l based on the pole structure in Eq. ( |72] |. This calculation 
involves a natural extension of the collective mode calcula- 
tions performed at the mean field level ll66[|67ll . but here one 
uses the fully self consistent jj* . In addition there has been 
work on the collective modes using NSR-2 which addresses 
an improved ground state which includes quantum fluctua- 
tions [47]. This, thus, goes beyond the mean field calculations 
of this earlier work, and quantifies the changes in the sound 
velocity. 

D. Alternative Schemes 

In this Review we have confined our attention to the two 
schools of BCS-BEC crossover which represent natural exten- 
sions of the seminal iflTifTsIl NSR and Leggett papers. There 



are alternate approaches which have been introduced into the 
literature. Most notable among these is a scheme associated 
with Zwerger ll90ll . Haussmann iFtsIi and their collaborators. 
The original work ifish could be viewed as a third alternative 
T-matrix scheme in which the pair propagator x{Q) appear- 
ing in Eq. ( |28T l involves two dressed Green's functions. In 
the context of work on high temperature superconductors, this 
scheme (and a closely related approach known as "fluctua- 
tion exchange" or FLEX) has been addressed by a number of 
different groups |[9ll l92l l93ll and there has been some contro- 
versy |91, 93, 9% about whether pseudogap effects naturally 
emerge. This approach has recently been extended [90] below 
Tc in somewhat the same spirit as the NSR-based schemes. 



V. DETAILED COMPARISONS 
A. Overview of Salient Qualitative Comparisons 

Because this paper is principally aimed at addressing theo- 
retical issues, we do not review the vast number of theory- 
experiment comparisons now in the literature. These are 
based on radio frequency spectroscopy, thermodynamics, col- 
lective modes and other techniques. Rather, here we address 
some of the major "milestone" issues which are often used 
to assess the general quality of a given BCS-BEC crossover 
theory. We begin with Table U which presents an overview 
of the two theoretical schools as summarized in Sections |III] 
and |IV] The first two lines characterize the behavior of the 
fermionic and bosonic dispersion as they appear in the re- 
spective "gap equations" of the two schools. As is consistent 
with the hierarchy of ground state wavefunctions in Section 
III Al one can infer that the NSR-based scheme approximates 
the fermionic contribution and focuses more directly on the 
bosonic contribution; it thereby arrives at a linear dispersion 
for the pairs. By contrast the BCS-Leggett school approxi- 
mates the bosonic contribution and focuses more directly on 
the fermionic dispersion, thereby incorporating pseudogap ef- 
fects into i?k- The order of the transition at Tc is second order 
in the BCS-Leggett scheme and first order lIsTll in NSR-based 
approaches. The finite T density profiles in a trapped gas will 
reflect this behavior and be rather smooth and featureless in 
the BCS-Leggett scheme ||95il while there will be derivative 
discontinuities and non-monotonic features |45] which reflect 
the condensate edge in the NSR based scheme. Similarly the 
first or second order of the transition will also show u p in 
the superfluid density within the BCS-Leggett 12ll [39l l95ll 
which displays smooth monotonic behavior or NSR based 
||48[ [ssl 89] scheme which shows multi-valed or discontinu- 
ous features at Tc- We point out that these spurious first order 
effects are also seen in the Bogoliubov theory for true bosons, 
as discussed in Appendix ICl 

Calculations of the critical velocity have been addressed 
within the BCS-Leggett school using Bogoliubov deGennes 
theory |55] and from Josephson junction studies |78]. Here an 
experimental comparison can also be made and the agreement 
II79II is reasonable. Table J] shows that in both schools the su- 
perfluid fraction in the ground state is 100 % in both schools. 
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NSR Based Scheme for general T 


BCS-Leggett Based Scheme for general T 


Fermionic Dispersion 
Rplnw T 


K=y(^k-M*)' + AL(r), 

3.pproxinistc treatment of fcrmions 


= V(ek - m)" + AL(r) + A2,(T) 


Bosonic Dispersion 
Rplnw T 


= c{T)q 


OPq — q^ /2M* , approximate treatment of 
bosons 


Order of Transition at Tc 


First order 


Second order 


Density Profiles at 
Unitarity 


Features indicating condensate edge 


Smooth, quasi-Thomas-Fermi 


Superfluid Density 


Multi-valued or Discontinuous at Tc 


Smooth and Monotonic at all T 


Calculations of Critical 
Velocity at T = 




Either From Vortex or Josephson Effect 

l^Q 1/^111 nf 1 r\nc 


T=0 Superfluid fraction at 
Unitarity 


100% 


100 % 


Order parameter 
collective modes 


LU = c{T)q 


UJ = c'(T)q 


Major Advantage of 
Ground State 


Captures physics of Bogoliubov theory, esp. 
good for BEC 


Allows spatial dependence. 

via Bogoliubov deGennes theory 



Table 1: Comparison of Conceptual Issues in the Two Different Theoretical Schools 



In the NSR IH IH and BCS-Leggett fs^] schemes the dis- 
persion of the order parameter displays the expected linear be- 
havior at long wavelengths. Finally we address the strengths 
of both ground states by noting that the NSR-based scheme 
captures the physics of Bogoliubov theory and should, thus 
be the quantitatively better ground state, particularly in the 
BEC limit. By contrast the BCS-Leggett scheme is the more 
flexible and allows a spatial dependence to be readily incorpo- 
rated in the form of Bogoliubov deGennes theory. Moreover, 
within the BEC, this BdG theory leads to a Gross Pitaevski 
picture of the ground state, which allows one to exploit a well 
established body of literature on true Bose systems. 



B. Comparison of Superfluid Transition Temperatures 

Figures[T]and|2]present comparisons of the superfluid tran- 
sition temperatures in the two schemes for the homogeneous 
situation and in a tra p. T he black lines correspond to the BCS- 
Leggett scheme il2lll3ll and the red lines are for the NSR ap- 
proach as obtained in Reference |45]. For the homogeneous 
case, it can be seen that there are only small quantitative differ- 
ences, while in the trapped situation the BCS-Leggett scheme 
leads to considerable lower Tc values slightly above unitarity. 
The root of the difference in the two calculational schemes 
lies physically in the fact that the BCS-Leggett scheme com- 
putes the transition temperature in the presence of a finite 
(pseudo)gap at Tc- In the NSR based scheme, these pair fluc- 
tuation effects do not appear as a pseudogap in the expression 
for Tc, but rather enter through corrections to the fermionic 
chemical potential /i* . 

Section ITlIBI presented simple arguments which show that 
the ideal gas asymptote for Tc is approached from below in the 



BCS-Leggett scheme, while it evidently is approached from 
above in the scheme of Nozieres and Schmitt-Rink. Both of 
these are mean field approaches and the behavior should not 
be compared with expectations |96] based on a critical fluctu- 
ation description of true Bose systems which clearly include 
other physical mechanisms. Indeed, the fact that at Tc there is 
a discontinuity in the NSR-based schemes suggests that this 
approach should be more suitable at T w 0, away from Tc- 

Because of the different approaches to the ideal gas asymp- 
tote, in a trap one sees from Figure |2] that the differences be- 
tween the two transition temperatures are more marked. The 
ideal gas asymptote is quickly reached in the NSR scheme 
very close to the point where l/fc^a « 1. In the BCS-Leggett 
scheme there is an extended regime at and on the BEC side of 
unitarity where Tc is rather constant, and the asymptote is only 
reached for 1/kpa considerably larger than its counterpart in 
the alternate school. 



C. Comparison Of Density Profiles 

Figure [3] presents a plot from Ref . | |46ll of the axial den- 
sity profiles in the BCS-Leggett ground state (dashed lines) 
as compared with the NSR-derived ground state (black lines) 
and the data points (shown in red) for ^Li. In axial profiles 
two of the three dimensions of the theoretical trap profiles 
were integrated out to obtain a one-dimensional representa- 
tion of the density distribution along the transverse direction: 
n{x) = Jdydzn{r)- Three different values of the magnetic 
field near unitarity are shown, and the upper and lower panels 
correspond to slight changes in the number of atoms, N which 
are assumed in the theoretical calculations. The figure shows 
that the agreement between theory and experiment is better 
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Figure 1: Comparison of Tc/Tp as a function of inverse scatter- 
ing length l/kpain a homogeneous system within the BCS-Leggett 
scheme L12i . il3i1 (black curve) and the Nozieres Schmitt-Rink |45] 
(red curve) scheme. The former has a maximum closer to unitarity 
and a dip close to the point where /j, changes sign. 




position (tim) 



Figure 3: Nozieres Schmitt-Rink based results for ground state pro- 
files. Shown is a comparison between experimental and theoretical 
axial density profiles from Ref. (4^ . Experimental data from Ref. 
(9^ (dots) are shown for three different values of the magentic field 
B tuning. Theoretical results for NSR theory 1 46] at T=0 (solid lines) 
and for BCS-Leggett theory (dashed lines) are shown for the same 
corresponding parameters. The upper (lower) panels refer to the es- 
timated number of atoms iV = 4 x 10^ (A'' = 2.3 x 10^). 




Figure 2: Comparison of Tc/Tf in the trapped case for the two 
schools, using the local density approximation. As in the previous 
figure, the red curve is for the NSR scheme |45] and the black curve 
for the BCS-Leggett approach H2i . il3i1 . 



for the smaller value of N. While the difference in the pro- 
files associated with the two ground states is not particularly 
dramatic, it should be stressed that this difference is reflected 
in rather large changes in the coefficient /3 discussed below. 
Overall the quantitative agreement between theory and exper- 
iment is seen to be better for the NSR-based ground state. 

In Figure|4]are shown density profiles at finite temperatures 
for the BCS-Leggett case, from Reference f95]. The experi- 
mental data and theory correspond to roughly T/Tp = 0.19. 
These profiles are estimated to be within the superfluid phase 
(Tc K, O.STp at unitarity). This figure presents Thomas-Fermi 
fits Wh to the experimental (0^) and theoretical d^') pro- 
files as well as their comparison (|4};), for a chosen Rtf — 
100 /iTO, which makes it possible to overlay the experimen- 
tal data (circles) and theoretical curve (line). Finally Fig. |4}l 
indicates the relative or root-mean-square (rms) deviations 
for these TF fits to theory. This figure was made in collabo- 
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Figure 4: (color online) BCS-Leggett results for temperature de- 
pendent profiles of (a) experimental one-dimensional spatial profiles 
(circles) and TF fit (line) from Ref. |97], (b) TF fits (line) to theory 
at T r; 0.7Tc ~ O.IQTf (circles) and (c) overlay of experimen- 
tal (circles) and theoretical (line) profiles, as well as (d) relative rms 
deviations (x^) associated with these fits to theory at unitarity. The 
circles in (b) are shown as the line in (c). The profiles have been 
normalized so that N = J n{x)dx — 1, and we set Rtf ~ 100 fim 



in order to overlay the two curves. 
T = 0.19Tf. 



X 



reaches a maximum around 



ration with the authors of Ref. ||97ll . To probe the deviations 
from a TF functional form, in Fig. |4}l, the (relative) rms de- 
viation, or x^, from the TF fits as a function of T is plotted. 
X^ increases rapidly below Tc and reaches a maximum around 
0.7Tc. Quite good agreement between theory and experiment 
is observed here in the finite T profiles. 
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D. Addressing (3 Effects 

One of the most widely used milestones for assessing 
crossover theories is the numerical value obtained for the co- 
efficient (3. At unitarity the chemical potential must scale with 
the Fermi energy with a coefficient of proportionality 



(93) 



In the BCS-Leggett ground state (3 ~ —0.41. By contrast 
experimental data IgSfl suggest an answer which is closer to 
Monte Carlo calculations |99(] (3 = -0.56. Calculations IH] 
based on NSR-1 yield j3 — —.545, which is in quite good 
agreement with experiment. In the NSR-2 scheme (based on 
the thermodynamical potential with the non-variational con- 
tribution included), the same good agreement with experiment 
(and Monte Carlo) was presented ll47i Issllsill . 

The BCS-Leggett- based scheme, as it has been imple- 
mented here, can be seen to ignore Hartree effects, just as 
is consistent with the ground state wavefunction. This is a 
shortcoming of the scheme and in the context of the formal- 
ism presented here, one can trace it to Eq. ( l34l i. This approxi- 
mation, in effect, includes only pairing correlations in the self 
energy. Correlations that are not associated with pairing, such 
as Hartree effects have been omitted. With the full T-matrix 
formalism as outlined in Section ITllAI it should be clear that 
this assumption can be avoided and is not fundamental to the 
physical picture presented here. However, dropping this sim- 
plification does lead to considerable numerical complexity. 
We note that the NSR-based theories both above and below 
Tc include these non-pairing correlations in a fairly automatic 
way. In some limited contexts, they have also been included 
in the BCS-Leggett based theory [30, 68, 100]. 

One can think of these omitted correlations as enter ing via 
Eq. ( |28]) when the pair susceptibility is assumed llOlll to in- 
clude only two bare Green's functions. Most recently, it has 
been shown that these "GoGo" correlations are responsible for 
some important physical observations in the context of gases 
which are so strongly polarized that superfluidity is driven 
away 1 102]. A bound state associated with the minority spins 



is found to occur I, l(j3[ 104ll in these highly imbalanced gases. 



whi ch is responsible lll05Lll06l] for anomalies in the RF spec- 
tra l|ip5]. 

In summary, it is possible to estimate the size of these 
Hartree corrections, if one goes beyond Eq. ( l34l ) and includes 
the effects deriving from "GqGo" correlations, noted above. 
This will have to be done in future calculations for better 
quantitative comparisons of various properties, including (3. 



E. Condensate Fraction and "Quantum Depletion" 

It is interesting to contemplate the concept of "quantum de- 
pletion" in a fermionic system, particularly as one approaches 
the BEC. It is generally believed that the BCS-Leggett the- 
ory is to be distinguished from that based on NSR (which is 
closer to Bogoliubov theory) because of the neglect of quan- 
tum depletion. However, because of the presence of unpaired 



fermions, the condensate fraction will automatically show less 
than 100% condensation, except in the deepest BEC. Simi- 
larly, in the BCS regime this condensate fraction is vanish- 
ingly small. Establishing the degree to which "quantum de- 
pletion" is present in a Fermi gas is, thus, a subtle issue. 

We begin with the BCS-Leggett ground state. Following 
earlier work fl08*], at T = 0, the pair wavefunction is defined 
as -Fk = (A^ — 2|c_kiCkT 1^)^ where is the fermion an- 
nihilation operator for a =t, i. It can be shown that in this 
ground state we have Fk = UkWk, where the coefficients are 
ul, vl = [l± (ek - ^i)/E^,]/2 and = V(ek - m)^ + A^. 
The condensate fraction at T = is 



^c = ^|Fk|2 = j dr\F{i 



(94) 



Here F{r) — ^^Fkexp(«k • r). This pair density reflects 
off-diagonal-long-range-order There have been a number of 
numerical calculations of this quantity over the entire BCS- 
BEC crossover |107] and the agreement with direct Monte 
Carlo schemes |[ 107.1 is not unreasonable, as will be summa- 
rized below. 

It is natural to try to extend this picture to finite tempera- 
ture, taking the quantity F]^ — T F{iuJn , k) = Wk^k [1 — 
2/(£^k)] as a measure of the pair density. We stress that this is 
not related to off-diagonal long range order, but rather con- 
tains the contributions from condensed and non-condensed 
pairs, through the decoupling of into A^^ and A^^. One 
has, thus. 



= A2[l-2/(£;k)]V4£^k, T^O. 



(95) 



To emphasize that there is no unique representation of the 
pair fraction away from the BEC limit, we note that Eq. ( l38b 
provides another natural decomposition We can rewrite this 
equation representing the total density of fermions n in the 
form 



n = 2ZA2 + 2^/(ek-M) 

k 

or equivalently 

n = 2ZAI + 2ZAl^ + 2 ^ /(ek - m), 



from which 



2ZAi 



2ZAi 



(96) 



(97) 



(98) 



can be obtained. There are three terms on the right hand side 
of Eq. ( |97] l. The second term corresponds to the density of 
fermions in the non-condensed pairs. The first term may be 
identified as Nc = 2ZA^^, representing an alternative way 
of quantifying the density of fermions in the condensate, and 
the third term may be identified as the density of remaining 
(unpaired) fermions, n/ — 2^^/(ek — /i). This decompo- 
sition is of interest, in part because it relates more directly to 
the decomposition of pairing contributions and free fermions 
introduced in the original NSR paper 
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NSR-1 


NSR-2 


BCS-Leggett 


"Answer" 


scatt. length ratio: asja 


2.0 


^0.55 


2.0 


0.6 (exact calc.) 




-0.545 


-0.59 


-0.41 


-0.55, (from experiment) 


atT=0, l/kpa = 2 




0.84 


0.99 


0.96 (Monte Carlo) 


atT=0, l/kpa = oo 




0.48 


0.69 


0.58 (Monte Carlo) 



Table II: Quantitative Comparisons among the different schools. References for each number are listed in the first row from left to right: 
Ref. jsl, Ref. (H], Ref. j?!, and Ref. |64]. For the second row from left to^right the references are: Ref. jH], Ref. (H, Ref. jH and 
Ref . f 98] . For the third row from left to right the references are Ref. Ref. 1^107] and Ref. 1 107]. Finally in the last row from left to right 
the references are Ref. 1,8^, Ref. 1107,1 and Ref. LI 07,1 . 



Recent calculations of Nc iiHlsol] have also been presented 
using the NSR-2 approach, where the fraction is found to be 
somewhat smaller than in the BCS-Leggett state. Importantly, 
the difference between these two results for Nc is viewed as 
a possible way to represent quantum depletion, which is nat- 
urally larger in NSR based theories as compared to the BCS- 
Leggett counterpart. 



F. Effects of First Order Transitions 

Essentially all NSR-based theories, as well as some which 
claim higher levels of consistency, report first order transitions 
lIsTll . These effects presumably originate in the same way as 
their counterparts in true Bose systems treated at the Bogoli- 
ubov Ii49i1 or Popov level. We outline the origin of these first 
order effects in Appendix |C] They lead to derivative disconti- 
nuities in the density profiles at the condensate edge lll09ll and 
are thus, not as problematic in the case of a trapped gases as 
compared to a homogeneous system. This is particularly the 
case in the BEC where bimodality is present and one would 
expect signatures of the condensate edge. 

Despite this theoretical framework, experiments show a be- 
havior which is far from first order. One of the most striking 
features about the unitary gases is that there is so little indica- 
tion of the phase transition and thus no evidence for first or- 
der behavior This is seen by noting the historical difficulties 
encountered in establishing whether a particular experiment is 
performed in the superfluid or normal phase. In the absence of 
population imbalance, the unitary gas profiles are featureless 
081] with no clear bi-modaUty or other indications of a conden- 
sate edge. Similarly, RF spectroscopic studies of the pairing 
gap show a smooth behavior [110] from high T to temper- 
atures well below Tc- As a consequence, Tc is difficult to 
identify, although important thermodynamical measurements, 
have indeed, indicated a phase transition [8]. 

These theoretically generated first order effects become 
even more difficult to reconcile with the fact that in BCS- 
BEC crossover, the pseudogap, which appears well above Tc 
leads to an even smoother transition than in strict BCS theory 
(which also is of second order). Thus a first order transition 



in systems undergoing BCS-BEC crossover can be viewed 
as somewhat problematic, except, perhaps if attention is re- 
stricted to a narrow temperature range. This points to an ad- 
vantage of the BCS-Leggett based approach where the density 
profiles are rather featureless and well fit to a Thomas Fermi 
form. A related advantage is that witho ut fir s t ord er transi- 
tions one can aiTive at a theoretical basis llllll llI2ll for adia- 
bat ic sw eep thermometry. This is an experimental technique 
10, ll 13il which has been rather widely discussed. Using the 
theoretically determined entropy, it is possible to arrive at rea- 
sonable estimates of a final temperature, based on an experi- 
mentally known initial temperature connected by an adiabatic 
sweep. 



G. Quantitative Comparisons 

Quantitatively, the NSR-based approaches appear to have 
some advantage, although there are variations depending on 
how the number equation is implemented (either via NSR-1 
or NSR-2). Table HH summarizes specific key numbers which 
have been used to assess the different schools. Listed in the 
first row is the scattering length ratio associated with the inter- 
boson and inter-fermion interactions, while the second row 
compares the quantity (3 defined in Eq. (|93l l which is associ- 
ated with the unitary limit. Finally, the third and fourth rows 
address the value of the condensate fraction in the ground state 
in the near-BEC and very deep BEC. This is relevant to quan- 
tifying the degree of quantum depletion. We have previously 
addressed our concerns about NSR-2 which, through Eq. dSST l. 
builds on inconsistencies associated with the fact that the gap 
equation does not satisfy the variational condition. Rather 
we argue in favor of the approach we call NSR-1 which uses 
Eq. ( l67l i. While there seems to be considerable interest in the 
community in comparing numbers such as those shown in Ta- 
bleini because of just these concerns about more fundamental 
issues, we are of the opinion that it may be premature to give 
too much weight to the numerical comparisons shown in Ta- 
ble Instead we attach greater importance to Table |T] for 
assessing a given theory and for indicating new directions for 
improvement. 
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VI. SUMMARY 

One of the major goals of this review has been to clarify 
the genesis of a large number of contributions to the theoreti- 
cal literature by associating them clearly with one or the other 
theoretical approaches to BCS-BEC crossover We stress that 
these two theoretical schemes are different in the ground state 
and in their thermal properties. One should, thus, avoid the 
tendency to present results from the Leggett BCS ground state 
and simultaneously use the Nozieres Schmitt-Rink calcula- 
tions for treating T 7^ aspects of BCS-BEC crossover. 

It was also our goal to summarize the major strengths and 
weaknesses within these two schools. While it clearly in- 
cludes bosonic degrees of freedom, there is a concern about 
the Leggett-BCS theory which concern derives from the fact 
that this approach does not yield a Bogoliubov-like or sound 
mode dispersion for the non-condensed pair excitations. For- 
mally, this is a consequence of the associated T-matrix de- 
scription of the (7 7^ pairs, which drops higher order terms 
and which are needed to couple the condensate and pair exci- 
tations. Although it has not been seen as yet 1 114], this sound- 
like excitation spectrum could show up in future experiments 
on unitary gases, particularly through power law dependences 
in thermodynamics. On the positive side, calculations in this 
BCS-Leggett phase are very tractable; one can readily handle 
inhomogeneities such as vortices (through the Bogoliubov de- 
Gennes approach); one can introduce trap effects, as well as 
population imbalance and address all temperatures T. 

For the NSR based schemes a comparably major problem is 
that there is no satisfactory mean field theory for the weakly 
interacting Bose gas which works at all temperatures. The 
NSR based school is based on this Bose gas mean field start- 
ing point and this introduces unphysical first order transitions 
which, at least around Tc, will interrupt the smooth crossover 
from BCS to BEC and limit the applicability of theory to spec- 
ified ranges of temperature. On the positive side, it is believed 
that this scheme, which works best at low temperatures, will 
produce the better ground state and allow more quantitative 
comparison with experiment at T w 0. 



With these weaknesses identified, each of these schools has 
a large agenda in hand for future research. In the short term 
the BCS-Leggett scheme should be readily extended to in- 
clude additional non-pairing contributions to the self energy 
(such as Hartree effects) which will make it more favorable 
for quantitative comparisons. Similarly, in the short term, the 
NSR-based approach should be extended to implement the in- 
clusion of Popov like correlations, and thereby include cor- 
rections to the BCS gap equation (which treats the fermions 
as non-interacting). In the longer term one would hope that 
NSR scheme (which approaches the crossover from the BEC 
end and oversimplifies the fermionic dispersion) and the BCS- 
Leggett scheme (which approaches the crossover form the 
BCS end and oversimplifies the bosonic dispersion), will ul- 
timately be unified. It is also to be expected that experiments 
will guide the way. 

Appendix A: DETAILS OF THE T-MATRIX CALCULATIONS 

The T-matrix is t-^{Q) = [C/^^ + x{Q)]- Near can it 
be expanded near Q = in the form 



tpg(q,rj) 



1 



Z{VL-q'^/2M* 



■iT. 



Q 



(Al) 



after analytic continuation (zil; ^ + iO+). The pair chem- 
ical potential fipair vanishes below Tc- The relaxation term 
iTq is neglected near T^. in most applications. Firstly we cal- 
culate Z. 



1 



0=0, 9=0 2A- 
The effective mass of pairs, M*, is given by 
1 1 d%Hq,0) 



(A2) 



2M* 



6Z 



q=0 



(A3) 



Here 



9=0 2A2 

^fiEl + ie^-p?)- 



}(Vkek)'] +4/'(ek-M)(Vkek)' 



(A4) 



Appendix B: DETAILS ON EQS. ^77) AND <80) 



We review how Eqs. ( ITTI i and (ISOl l are derived following 
iQ. The Dyson's equation, Eq.(32) in Ref.iH, is 
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Gl2{K) = Gi2(i^)=G?i(i^)[Sll(i^)Gi2(i^)+Si2(i^)G22W]. 

I 



(Bl) 



Following the approximation shown in Eq. dTSl l. the self- 
energy becomes 



A2 

pg 



(B2) 



In the BEC limit, Asc/Im*! ^ 1 ^i^^ we assume this also 
holds for Apg/|^*|. Then Gii{K) 

A 2 

itOn + ^k - 



- Ck - f + Ck - ) - A? 



2 

k ' "pg- 



(zc^„ + ek)KHek + A^„) 



ii^n + 'Ck 



^2 +^2^A2g+A 

^^ri + Ck 
^2 +£2 + A2 + A2 ■ 



(B3) 



Here we made the approximation, which is valid in the BEC 
limit, 



+ Ck + Ha 



= 1 - 



P3 



w2+£2+A2 



1. (B4) 



^k ' "pg 

The off-diagonal fermion Green's function G12 (K) 
A^^ 

A,,(c^2+£2) 



(t^ue^ + A2^)2 + ALK + a 

A^e 

,,,2 1 C2 I o X 2 I A 2 J _ 

'^n + Sik + ^^pg + ^sc + (i^2_|_^^) 

A,, 



a;2+£2+2A2^ + AL 



(B5) 



Note that the corresponding gap equation, Eq. JsTT l. derived 
from this expression is different from the BCS gap equation 
when is defined as ^£2 + /^2^ + ^2^_ 



Appendix C: FIRST ORDER TRANSITIONS IN BOSON 
MEAN FIELD THEORIES 



It is well known 11501 ll 1511 that mean field theories of the 
weakly interacting Bose gas are associated with unphysical 



first order transitions. It is, thus, often argued that these the- 
ories should only be applied at temperatures much lower than 
Tc- Since the same issues arise with BCS-BEC crossover the- 
ories of the extended Nozieres Schmitt-Rink school (based on 
the Hartree-Fock-Bogoliubov or Popov approximations), it is 
useful to understand the physical origin of these first order ef- 
fects. 

We summarize here the central issues which lead to first 
order transitions: 



1 . The BEC transition temperature predicted by mean field 
theories is the same as the BEC temperature of an ideal 
Bose gas, Tg^c- 



2. Below Tg^p, interaction effects are found to suppress 
thermal excitations. This suppression arises from inter- 
action effects in the dispersion relation which lead to a 
systematic increase in the excitation energy, relative to 
the non-interacting gas. In addition there is a change in 
the phase space weighting factor In combination, these 
two effects importantly yield a smaller fraction of non- 
condensed bosons (or a larger condensate fraction). 



3. As a consequence, if one plots the condensate fraction 
obtained from generic mean field theory as a function 
of T, one sees that it tends towards Tg^^ by overshoot- 
ing and then bending back towards Tg^^ at the highest 
temperatures below the transition. This double valued 
behavior is then associated with a first order transition. 



Fig Hshows the condensate fraction (solid line) as a func- 
tion of temperature as obtained from the Popov approxima- 
tion. The bend-over which indicates a first order transition 
can be seen clearly. As a comparison, the condensate fraction 
of a non-interacting gas of bosons is also presented (dashed 
curve). Here one sees a smooth second order phase transition. 

We can take these ideas over to the BCS-Leggett approach 
to BCS-BEC crossover (which is the only case where a first 
order transition is not seen). It is rather straightforward to 
see the analogies with the Bose gas through the gap equation: 
A2 = A2^ + Apg in conjunction with Eq. ( [35] l (and in some 
situations also with Eq. (IJTTi). Here, too, following the inter- 
acting Bose gas logic we will also end up with an unphysical 
first order transition which means that the smooth crossover at 
finite T is interrupted for some range of temperatures below 
and near Tc. 
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Figure 5: Condensate fraction no/*i as a function of dimension- 
less temperature T /T%^q. Solid and dashed lines correspond to 
results obtained from the Popov approximation and the theory of 
non-interacting bosons, respectively. 




Here Gr = [G(7^(if ) - T.{K)]-^ is the full fermion Green's 
function. The full vertex is obtained by choosing a set of di- 
agramms consistent with gauge invariance. These are shown 
in Fig. |6] and specified in the caption. Now we would like to 
show that this set of diagrams satisfy the Ward identity and 
does not contribute to the Meissner effect if the correct form 
of number equation is used. 

The expressions of those diagrams in Fig.|6]are 



\{K,K+) = (2k + q,l), (D2) 

MT{K,K+) = Y.*{P)Go{P - K+) X 
p 

X{P^ K+,P~ K)Go{P~ K), 
AL{K,K+) - -2j2tiPHP+)Go{P - K) X 

L,P 

Go{P - L)Go{L)X{L, L+)G„{L+). 



Figure 6: The full vertex is approximated by these diagrams. The 
first one on the right hand side is the bare vertex, the second one is 
the "MT" diagram, and the last one is the "AL" diagram. Hollow and 
solid dots denote full and bare vertices. Solid lines and wavy lines 
correspond to propagator of non-interacting fermions and t-matrix, 
respectively. 



Appendix D: WARD IDENTITY ANALYSIS OF THE 
NUMBER EQUATION 

The Ward identity is given by 



Q ■ XiK, 
Q ■ A{K, 



) = G,\K) - G^\K+), 



G-\K) - G-\K+). 



(Dl) 



The factor 2 in the AL diagram comes from the fact that the 
vertex can be inserted in one of the two particle propagators 
in the T-matrix. By taking inner produc with Q they becomes 



Q-MT{K,K+) = J2tiP)Go{P - K+)[G,\P - K+) 



G^\P - K)]Go{P - K) 



(D3) 



Q ■ AL{K,K+) = -2Y,t{P)t{P+)Go{P - K)Go{P - L)Go{L)[G^\L) - Go\L+)]Go{L+), 

= -2 ^ t(P)i(P+)Go(P - K)[x{P+) - x{P)]. 
p 

= 2[S(i^+) - ^{K)]. 



In deriving these results, Eqs. (ID2l) and the identity x{P+) ^ 
X{P) = t-'^{P+) - t-\P) are useful. Therefore Q ■ {MT + 
AL) = T,{K+) - ^{K). It is straighforward to show that the 
approximated vertex A = A + MT + AL satisfies the Ward 
identity. Since these diagrams are normal state diagrams, one 



I 

can take the limit Q ^ and obtain 

[MT{K, K) + AL{K, K)]^, 



(D4) 



(D5) 



This is an important identity in the derivation of the absence 
of Meissner effect in the normal state. 

The Meissner effect occurs if the static response kernel does 
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not vanish, which is equivalent to the existence of a superfluid 
density. To show that the approximation for the full vertex 
does not contribute to the Meissner effect, it suffices to show 
that in the normal state 



(-) - Pxx(O) = 0. 
Here the density is calculated as (a, /? = x, z) 



(-) = 2E 

\mJa8 ^ 



^771 y cc8 



-G{K) 
d^k dG{K) 



(D6) 



(D7) 



K 



K 



We assume that surface terms can be neglected. The current- 



current correlation function at Q = is 



PaM = -2Y,G\K)[X{K,K) + MT{K,K) + 



K 



AL{K,K)]0\{K,K)„ 



K 



dK„ 



dK, 







,(D8) 



where we used X{K, K)^ = 4^ 



9K„ 

Thus the two contributions cancel each other and there is 
no Meissner effect in the normal state. 
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